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Abstract 

The cross sections for the reaction e + e~ — > I I q q and for similar four fermion production processes 
at LEP 1, LEP 2 and the NLC are calculated. Due to the extraordinary symmetry properties 
of the process, very compact analytical formulae describe the double differential distributions in 
the invariant masses of the 11 and qq pairs. The total cross sections may be obtained with two 
numerical integrations. 



email: bardindy@cernvm.cern.ch, leike@cernvm.cern.ch, riemann@ifh.de 

" Supported by the German Federal Ministry for Research and Technology under contract No. 05 GMU93P. 



1 Introduction 



At LEP 1 energies, few events of electron positron annihilation into four fermions, e + e~ — > Af , 
have been observed jl] although 4/ production is extremely suppressed compared to fermion pair 
production. The 4/ production is also interesting as a background contribution to light Higgs 
searches with the associated ZH production. The strong kinematical suppressions are hoped 
to become overcompensated by the rising number of Z bosons being produced at LEP 1. 

Above the production thresholds for gauge boson pairs, four fermion production is predicted 
by the Standard Model to be one of the most frequent annihilation processes. It will allow to 
study further details of the gauge boson properties; and the search for Higgs boson signals from 
the Bjorken process will be one of the challenges of LEP 2. 

Several Monte Carlo approaches to the complete description of 4/ production have been de- 
veloped ||, |^, |6| while we are performing a program of semi- analytical calculations [0, [|, 0. 
Although the treatment of even the simplest final states is technically involved we got encouraging 
results for off shell W pair production |]7], || and the Bjorken process (off shell ZH production) ||. 
Some general remarks on 4/ production and numerical comparisons may be found in || |10"1| . 
The neutral current (NC) channels for 4/ final states are classified in table 1. The four different 
event classes are discussed in ||. Higgs boson exchange has been excluded; this is a very good 
approximation if one doesn't perform a dedicated Higgs search ||. 

In this article, we deal with the simplest event class marked in boldface. It comprises fi- 
nal states which do neither contain electrons (positrons), electron (anti)neutrinos, nor identical 
fermions and cannot be produced by charged current interactions. 

Our calculation covers the following observable final states: 

(i) [fifi,Tf], 

(ii) [11, bb], [11, cc), [11, {uu + dd + ss + cc + bb)), 

(iii) [bb,cc], 

where I = fi,r. Other reactions of table 1, e.g. those with u T and pairs, may also be calculated 
with the below formulae but are not observable due to the admixture of other final states, here 
e.g. v e pairs. The latter belong to another, much more complicated type of reaction. 

In section 2, the calculation is shortly described. Section 3 contains the main result: the 
analytical formulae for the invariant fermion pair mass distributions. In section 4 we discuss 
numerical results. 
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Table 1: Number of Feynman diagrams for 'NC type final states. 
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2 Feynman diagrams, phase space, and cross section 



The reaction 



e e 



ll + bb 



is representative for the class of reactions considered here. It is described by 24 Feynman diagrams 
with six different topologies which we will call crab diagrams and deer diagrams; see figures 1 
and 2. If there are two quark pairs in the final state, the deers may also contain gluon exchange. 
This will be discussed later. 

We parametrize the eightdimensional four particle phase space as follows: 



dT 



i=l ^Pi i=l 



^\{s,s llS2 ) ^\( Sl ,mi,m 2 2 ) yJX(s 2 ,mi,ml) 
= 2ir- as ids 2d cos Uailidi 2 2 , 



8s 8si 
with the usual definition of the A function, 



8s 2 



(2) 



X(a,b,c) = a 2 + b 2 + c 2 - 2ab - 2ac - 26c, 
A = X(s,si,s 2 ). 



(3) 
(4) 







Figure 1: The crab diagrams. 





Figure 2: The b-deer diagrams. The 1-deers may be obtained by interchanging the leptons with 
the quarks. 
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In (|2|), the rotation angle around the beam axis has been integrated over already. Variables k\ 
and hi are the four-momenta of electron and positron and Pi,P2,P3,P4 are those of the final state 
particles fi,fi,f2, f2 with pf = mf. All fermion masses are neglected compared to the invariants 
s, Si, and s 2 Q: 

s = (ki + k 2 ) 2 , si = (pi + p 2 ) 2 , s 2 = (p 3 + Pi) 2 - (5) 

The angle 9 is located between the vectors (pi + P2) and k\. The spherical angle of p\ [p^) in the 
rest frame of the compound [fifi\ ([/2/2D is ^i (^2): dVLi = d cos Oidfa. The kinematical ranges 
of the integration variables are: 

(mi + m 2 ) 2 < si < (a/s - m 3 - m 4 ) 2 , (m 3 + m 4 ) 2 < s 2 < - v^i) 2 ' 

— 1 < COS0, COS^i, COS 6 2 < 1 , < 2 ,0i < 27T. (6) 

We are interested in analytical formulae for distributions in invariant masses of fermion pairs. 
Thus, we have to integrate analytically over the five angular variables in (0). The matrix elements 
squared have been determined by two independent calculations. One of them used the symbolic 
manipulation program FORM ||11|| . The other one determined the squared matrix elements with 
CompHEP JT^I and continued then with the angular integrations in a different FORM program. 
The total cross section may be written as a sum of six contributions: 



- I- dS2 ^ a~7^~ 

In (^) we allow for lower cuts on the invariant masses. Index k = 1 corresponds to the square of 
the sum of the crabs, and k — 2, 3 to that of the b-deers and 1-deers. Further, k = 4, 5 and 6 
correspond to the interferences between 1-deers and b-deers, between crabs and b-deers, and 
between crabs and 1-deers, respectively. 



3 The cross section 

All the diagrams which mediate process ([!]) have a common topology. They contain three fermion 
antifermion pairs connected by neutral gauge bosons. Two of the fermion pairs couple to a gauge 
boson while the third one couples to two gauge bosons with a fermion propagator between them. 
In this sense, there is only one generic topology, the deer topology, and the crab diagrams may 
be considered as the e-deers in a crossed channel. 

Our choice of coordinates takes the above symmetry into account. All the 300 interferences 
between the 24 diagrams may be reduced to two generic ones. This may be achieved by grouping 
the diagrams in the calculation as is indicated in the figures, i.e. in sums of deers of the same 
type. A tremendous reduction of terms arises during the calculation leading to an extremely 
compact result. 

The first generic interference is known from off shell Z boson pair production ||. It is the square 
of the crab diagrams or generally of same-type deers. The three cases differ from each other by 
permutations of the invariant masses and couplings. The second generic kinematic function is due 
to the interference of different types of deers. Again, three such cases are related by permutations. 

1 An exception arises if there are photon propagators 1/s; involved and if no cut is applied. In this case, the 
finite masses of the fermions yield non- logarithmic contributions of order 0(1). These may be taken into account 

by replacing in (|^) the « /l — 4m^/si by \h-~ 4m^/sj (1 + 2mj/si) (instead of neglecting rrij/si). 
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Hence, the complete cross section can be written as a sum of six basic interferences between the 
three basic sets of diagrams. 

The cross section from the two crab diagrams alone reads as follows: 

d 2 a 1 (s;s 1 ,s 2 ) V\ . , , . . 
1—, = — ^ 422 (e, s;Z,Si;6,s 2 )&422(s;si,s 2 ) • (8) 

dSidS2 7TS Z 

This expression factorizes into a function of couplings and a kinematical ^-function and is sym- 
metric in Si and s 2 ; the latter depends exclusively on the virtualitiesQ: 

S 2 ~t~ f.S'i — ^2)^ 

^422(5; si, s 2 ) = — £(s;sx,s 2 ) -2, (9) 

S — Si — s 2 

with 

rf \ 1 , s-si-s 2 + v / A 

£(s; si, s 2 ) = -7= m =. (10) 

V A s — S\ — s 2 — v A 

The couplings and gauge boson propagators are contained in the following C-factor: 

C*422(e, s; fi, Si] f 2 , s 2 ) = 9fte V] — — p — - — — p — - „ , — 77777 — \ 

( 67r ) v^,l£^= 7 ,z D ^ (si) D ^'( s 2) £>y fe 0i) £>y ; (s 2 ) 

x [L(e, V^)L(e, ^)L(e, V})L(e, VJ) + i2(e, K<)i2(e, K fc )i2(e, ^)it!( e , V x )\ 

x [L(/i, V fc ) + Rifr, Vi)R(f u V k )\ N c {h) 

x [L{f 2 , V 3 )L(f 2 , Vi) + R{f 2 , V 3 )R(f 2 , Vi)] N c (f 2 ). (11) 

The conventions for the left- and right-handed couplings between vector bosons and fermion / 
are: 

L(/, 7 ) = B(/,7) = ! f £ , (12) 

LiS - Z) = 4^K- 2 ^)' fl ( / ' Z ) = 4^(- 20 ^)- < 13 > 

Here, we use e = V^ra, Q e = —1, J| = — |, and s 2 ^ = 1 — M^/Mf . Further, we choose a/7r 
= v^M^s^G^/tt 2 . The numerical input for the figures is G M = 1.16639 x 1(T 5 , M w = 80.220 
GeV, M z = 91.173 GeV, T z = 2.497 GeV. The colour factor N c (f) is equal to unity for leptons 
and three for quarks. The propagators are 

D v (s) = s-Arfr +Vs"r v (s). (14) 

For the photon, it is M 7 = T 7 = 0. The Z width function depends on the decay channels which are 
open at a given energy; at the Z peak and above, a good approximation is T z (s) = y/sTz/Mz- 
The D* is the complex conjugate of the propagator D. The resonating behaviour of the cross 
section depends on the arguments of the function C4 22 . With e.g. all Vi = Z in (|TT|), one selects 
the off shell Z pair production diagrams squared (see figure 1), whose contribution is proportional 
to two Breit Wigner resonance factors, 



C422 = • • • + 2 L (e,Z) + R e, Z) - , 15 

J 7T L» Z (Si) 2 7T .L>£(S 2 J 



2 This function may be found in the article p3| on electromagnetic pair production and also in |L4| where extra 
neutral gauge boson production is studied. 
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with 



l im I^L = 5 ( s _m 2 ). 



(16) 



r->o tt\D(s)\ 2 

The cross section contributions a 2 from the square of the b-deers and <7 3 of the 1-deers are: 



d 2 (T 2 {s;s 1 ,s 2 ) VX . i \r i \ 

■ = — 5-0422(0, s 2 ;e,s;/,si)C/ 4 22(s2;s, si), 



ds\ds 2 



7TS Z 



d 2 a 3 (s;s 1 ,s 2 ) VX 

= — ~C 422 (Z, si; b, s 2 ; e, s)£ 422 (si; s 2 , s). 



ds\ds 2 



7TS Z 



(17) 
(18) 



In these contributions the potential resonance behaviour is in the variables s and s± or s 2 , respec- 
tively. 

Besides the above moduli squares, there are interferences among the three different groups 
of diagrams. The interference between 1-deers and b-deers is again symmetric in the last two 
arguments: 



d 2 a 4 (s; Sl ,s 2 ) y/X 

— ^233l e ? S, t, Si, 0, s 2 )y 233 {s, Si, S 2 J, 



ds\ds 2 



7TS Z 



(19) 



with the kinematical function £233(5; Si, s 2 ), 

£233(5; si,s 2 ) = ^|>C(s 2 ;s,Si)£(si;s 2 ,s)4s ssi(s - Si) 2 + ss 2 (s - s 2 ) 2 + Sis 2 (si - s 2 ) 2 



+ (s + Si + s 2 



£(s 2 ; s, si)2s (s - s 2 ) 2 + si(s - 2si + s 2 ) 
+ £(si; s 2 , s)2s (s - si) 2 + s 2 (s + si — 2s 2 ) 
+ 5s 2 - 4s(si + s 2 ) - (si - s 2 ) 2 



(20) 



For A — * it remains finite, £233(5; Si, s 2 ) — > —3 (3s — s 1 — s 2 )/ [(s + Si — s 2 )(s + s 2 — si)]. The 
C factor with couplings and propagators is: 



C 23 3(e,s;/i,si;/ 2 ,s 2 ) = 3?' 



E 



(6vr 2 ) 2 y^,tW ^(5) ^(5 2 ) L>^(s) ^(Si) 
x [L(e,V^L(e,V k ) + R{e,V^R(e,V k )] 

x [L(/i, VJ)L(/i, Vj)L(h, Vt) - VJ)i2(/i, Vj)R(fi, V t )\ N c (fi) 
x [L(/ 2 , Vj)L{f 2 , V k )L(f 2 , Vt) - R(f 2 , Vj)R(f 2 , V k )R(f 2 , V{)] N c (f 2 ). 

(21) 

Both final fermion traces must couple to an axial current. Otherwise the contribution vanishes 
due to the Furry theorem. 

The remaining two interferences <r 5 among crabs and b-deers and cr 6 among crabs and 
1-deers are: 



d 2 a 5 {s; Sl ,s 2 ) VX ( , , , 

-j— i = — pC , 23 3(/,si;o,s 2 ;e,s)£ 2 3 3 (si;s 2 ,s), 

ds\ds 2 ns 2 



(22) 



d 2 a 6 (s;s 1 ,s 2 ) VX 



ds\ds 2 



7TS Z 



C 233 (o, s 2 ; e, s; /, Si)£ 233 (s 2 ; S, Si) 



(23) 
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The interferences of different-type deers may become resonating only in the first argument of 
the function C 233 (fi, si; / 2 , s 2 ; / 3 , s 3 ). 

At the end of this section, we would like to comment on processes with four quarks in the 
final state. In this case, besides photons and Z bosons also gluons may be exchanged in the deer 
diagrams while the crab contribution o~\ remains unchanged. The gluonic contributions to 05 and 
<j 6 vanish due to the colour trace, but a 2 , 03, and 04 get additional contributions. The sum over 
gauge bosons in the C-functions in QT7|) — (|I5D extends now also over gluons. The gluon couplings 
to quarks are 

L(q, 9) = R{q, 9) = lv^~ s , (24) 

and the gluon propagator looks like the photon propagator. Further, the colour factor N c (fi)N c (f2) 
has to be replaced by a factor of (AT? — l)/4 = 2 for all interferences with two gluons and by a 
zero if only one gauge boson is a gluon. 

If one observes two hadron jets without flavour tagging, then there are also contributions from 
diagrams with a three gluon vertex which have to be added incoherently. If the quark types may 
be tagged, such contributions are of higher order and may be neglected. 



4 Results 

Numerical comparisons have been performed for e + e~ — > fi + fi~bb and e + e~ — > Vpppbb at LEP 2 
energies with a Monte Carlo approach M; see also (TO). The agreement is within 0.2% and is 



limited by the numerical accuracy of the Monte Carlo program. 

There are two different energy regions of physical interest. At LEP 1, among others the 
cases (i) and (ii — 3) (see Introduction) are observed and case (ii — 1) is searched for as a Higgs 
boson signal HQ. While, at higher energies the most interesting channel is ( i i — 1 ) . 

In figure |3], the cross sections of cases (i), (ii — 1) , and (ii — 3) are shown at LEP 1 energies. 
We assume cuts being applied to the final state fermion pairs. For illustrational purposes, cross 
section (i) is shown also without cut. The cross section is sensitive to this due to the diagrams 
with photon exchange which are extremely enhanced when their virtualities approach to Am'j. 
The cross sections are peaking at ^/s = Mz, if the non-resonating gauge boson is a photon with 
a small virtuality. In partial widths of the Z boson into 4/ final states are systematically 
studied. Having in mind that the production cross sections are proportional to them at the Z 
peak, one may compare corresponding ratios with each other. We have done this in the massless 
limit for the cases [up,, rf), [pp, bb], [pp, (uu + dd+ ss + cc + bb)] and found agreement within 1 % 
which is the accuracy of the numbers quoted in [[|. 

Figure £| shows final states (i), (ii — 1) , and (iii — 1) in the energy region y/s ~ 100 — 500 GeV. 
Again, in one case the influence of the cuts on the invariant masses is exhibited. As to be expected, 
there are pronounced thresholds at the onset of on shell Z pair production. At higher energies, 
the cross sections fall monotonically. In the figure, we show also a four quark final state. The 
strong coupling constant a s is a function of the virtualities si or s 2 (and in some interferences of 
both) and may not be reasonably chosen here. For the illustrational purposes, we fixed it somehow 
arbitrarily at a s = 0.2; this corresponds to smaller values of the gluon virtualities which give the 
dominant contributions. 

We also should mention that there are substantial radiative corrections to 4/ production. The 
bulk of them is known to arise from initial state photonic bremsstrahlung. It may be easily taken 

3 We did not perform a dedicated comparison with the experimental results at LEP 1 jjj since the applied cuts 
are different from ours; our predictions agree roughly with them. 
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Figure 3: Total cross sections a(e + e — > Af) at LEP 1 energies. The following cuts are applied: 
Egg, E„n > 2 GeV, E Tf > 2m T , E cZ > 5 GeV, E b i > 20 GeV. 



Figure 4: Total cross sections a(e + e — > 4/) at £ and 7VLC energies. The following cuts 
are applied: E^ > 2 GeV, E Tf > 2m T , E cS > 5 GeV, E h i > 20 GeV. 



into account with a convolution formula which is well known from the Z line shape [||, |]. These 
corrections amount to - (9(15%) around the peaks and are smaller elsewhere, becoming positive 
(and cut dependent) in the radiative tail regions. Since all these effects are well known we do not 
go into numerical details on them. Our Fortran program 4f AN may take them into account. 

To summarize, we performed the first complete semi-analytical calculation of neutral current 
four fermion production. We obtained, for the simplest topology, extremely compact analytical 
results for the two-dimensional invariant mass distributions. The remaining two integrations are 
fast and numerically stable. Polarized beams and the inclusion of heavy neutral gauge bosons are 
also described by our formulae. 
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